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Abstract 

We propose a fairly simple and natural extension of Stollmann's lemma to correlated 
random variables, described earlier in [5|. This extension allows (just as the original Stoll- 
mann's lemma does) to obtain Wegner-type estimates even in some problems of spectral 
analysis of random operators where the Wegner's lemma is inapplicable (e.g. for multi- 
particle Hamiltonians) . 

To the best of author's knowledge, such an extension seems to be original. However, the 
author will appreciate any reference to articles or preprints where similar results are proved. 

1 Introduction 

The regularity problem for the limiting distribution of eigen-values of infinite dimensional 
self-adjoint operators appears in many problems of mathematical physics. Specifically, con- 
sider a lattice Schrodinger operator (LSO, for short) H : £ 2 (Z d ) — > f 2 (Z d ) given by 

(H^)(x) = Tp(v) + V(x)iP(x); x,yeZ d . 

y: \y-x\ = l 

For each finite subset A C Z d , let Ef,j = 1, . . . , |A|, be eigen-values of H with Dirichlet b.c. 
in A. Consider the family of finite sets Al = [-L, L] d (~)Z d and define the following quantity 
(which does not necessarily exist for an arbitrary LSO): 

If the above limit exists, k(E) is called the limiting distribution function (LDF) of e.v. of 
H. One can easily construct various examples of the function V : 1 d — > R (called potential 
of the operator H) for which the LDF does not exist. One can prove the existence of LDF 
for periodic potentials V , but even in this, relatively simple situation existence of k(E) is 
not a trivial fact. 

However, one can prove existence of k(E) in a large class of ergodic random potentials. 
Namely, consider an ergodic dynamical system (Q, T, P, {T x , x £ 1 d \) with discrete time 
1 d and a mesurable function (sometimes called a hull) v : Q — * R. Then we can introduce a 
family of sample potentials 

V(x,U!) = v(T x lj), x £ Z d , 
labeled by w 6 fi. Under the assumption of ergodicity of {T x }, the quantity 
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is well-defined P-a.s. Moreover, k(E, ui) is P-a.s. independent of u), so its value taken for a.e. 
u> is natural to take as k(E). In such a context, k(E) is usually called integrated density of 
states (IDS, for short). It admits an equivalent definition: 



where / G e 2 (Z d ) is any vector of unit norm, and II( — oo, E](H(lu) is the spectral projection 
of H(u>) on (— oo, E\. The reader can find a detailed discussion of the existence problem of 
IDS in excellent monographs by Carmona and Lacroix t 4! and by Pastur and Figotin |16| . 

It is not difficult to see that k(E) can be considered as the distribution function of a 
normalized measure, i.e. probability measure, on R. If this measure dK(E), called measure 
of states, is absolutely continuous with respect to Lebesgue measure dE, its density (or 
Radon-Nikodim derivative) dK(E) / dE is called the density of states (DoS). In physical lit- 
erature, it is customary to neglect the problem of existence of such density, for if dK(E) / dE 
is not a function, then "it is simply a generalized function". However, the real problem is 
not terminological. The actual, explicit estimates of the probabilities of the form 



for LSO H\ L in a finite cube Al of size L, for small e, often depend essentially upon the 
existence and the regularity properties of the DoS dk(E) /dE. 

Apparently, the first fairly general result relative to existence and boundedness of the 
DoS is due to Wegner [TH] . Traditionally referred to as Wegner's lemma, it certainly deserves 
to be called theorem. 

Theorem 1 (Wegner) Assume that {V(x,lo),x £ Z d } are i.i.d. r.v. with bounded density 
pv(u) of their common probability distribution: \\pv\\oc = C < oo. Then the DoS dk(E)/dE 
exists and is bounded by the same constant C . 

The proof can be found, for example, in the monograph [4]. 

This estimate and some of its generalizations have been used in the multi-scale analysis 
(MSA) developed in the works by Frohlich and Spencer [T3], Frohlich, Spencer, Martinelli 
ans Scoppola [14] . von Dreifus and Klein [llj . |12j . Aizenman and Molchanov p], and in a 
number of more recent works where the so-called Anderson Localization phenomenon has 
been observed. Namely, it has been proven that all e.f. of random lattice Schrodingcr 
operators decay exponentially at infinity with probability one (for P-a.e. sample of random 
potential V(ui)). Von Dreifus and Klein 12 proved an analog of Wegner estimate and used it 
in their proof of localization for Gaussian and some other correlated (but non-deterministic) 
potentials. The author of these lines recently proved, in a joint work with Yu. Suhov [9], 
an analog of Wegner estimate for a system of two or more interacting quantum particles 
on the lattice under the assumption of analyticity of the probability density pv(u), using a 
rigorous path integral formula by Molchanov (see a detailed discussion of this formula in the 
monograph 4 ). In order to relax the analyticity assumption in a multi-particle context, V.C. 
and Yu. Suhov later used ([10]) a more general and flexible result guaranteeing existence 
and boundedness of the DoS: the Stollmann's lemma, which we discuss below. 

In the present work, we propose a fairly simple and natural extension of Stollmann's 
lemma to correlated, but still non-deterministic random fields generating random potentials. 
To the best of author's knowledge, such an extension seems to be original, although very 
simple. However, the author will appreciate any reference to published papers or preprints 
where the same or similar result was mentioned and proved. Our main motivation here is to 
lay out a way to interesting applications to localization problems for multi-particle systems. 



k(E) = E[(f,Il ( - oo,E](H(u;)f)], 
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2 Stollmann's lemma for product measures 

Recall the Stollmann's lemma and its proof for independent r.v. Let m > 1 be a positive 
integer, and J an abstract finite set with |J|(= cardj) = m. Consider the Euclidean space 
R J = R m with standard basis (ei, . . . , e m ), and its positive quadrant 

R+ = | q G R J : q 3 > 0, j = 1, 2, . . . , m } . 

For any measure /i on R, we will denote by fi m the product measure pi x ■ ■ ■ x fi on R J . 
Furthermore, for any probability measure jx and for any e > 0, define the following quantity: 

s(p, e) = sup / dju.(t) 

a 

and assume that s(fi, e) is finite. Furthermore, let /i m_1 be the marginal probability distri- 
bution induced by p m on g' = (q2, . . . , g m ). 

Definition 1 Let J be a finite set with \J\ = m. Consider a function $ : R J — » R on R J 
which we will identify with R m . Function $ is called J -monotonia if it satisfies the following 
conditions: 

(f) for any rel+ and any q G R m ; 

$(g + r) >$( g ) ; (1) 

(2) moreover, for e — ei + ■ ■ ■ + e m G R m , /or an?/ g G R m and /or any t > 

$(g + i-e)-$(g) >l. (2) 

It is convenient to introduce the notion of J-monotonic operators considered as quadratic 
forms. In the following definition, we use the same notations as above. 

Definition 2 Let ri be a Hilbert space. A family of self-adjoint operators B : TC x R J — » Ti 
is called J-monotonic if, for any vector f £ Ti. with \\f\\ = 1, the function "!>/ : R J — > R 
defined by 

*f(q) = (B(q)f,f) 

is monotonic. 

In other words, the quadratic form Qs( 9 )(/) := (B(q)f,f) as function of q G R J is 
non-decreasing in any qj, j — 1, . . . ,\J\, and 

(B(g + I - e )/,/)-(B(g)/,/)> I .||/|| 2 . 



Remark 1 By virtue of the variational principle for self-adjoint operators, if an operator 
family H(q) in a finite- dimensional Hilbert space Ti is J-monotonic, then each eigen-value 
°f B(q) is a J-monotonic function. 

Remark 2 If H{q), q G R J , is a J-monotonic operator family in Hilbert space Ti, and 
Ho : Ti — > Ti is an arbitrary self-adjoint operator, then the family Ho + H(q) is also J- 
monotonic. 
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This explains why the notion of monotonicity is relevant to spectral theory of random 
operators. Note also, that this property can be easily extended to physically interesting 
examples where Ti has infinite dimension, but H(q) have, e.g., compact resolvent, as in the 
case of Schrodinger operators in a finite cube with Dirichlet b.c. and with bounded potential, 
so the respective spectrum is pure point, and even discrete. 

Theorem 2 (Stollmann, |17j ) Let J be a finite index set, \J\ — m, /i be a probability 
measure on R, and /i m be the product measure on R J with marginal measures fj,. If the 
function <& : R J — > R is J-monotonic, then for any open interval / Clue have 

V m {q: <f>(q)el}<m-s(p,\l\). 

Proof. Let I — (a,b), b — a = e > 0, and consider the set 

A = {q: $(?)<«}. 

Furthermore, define recursively sets Aj, j = 0, . . . , m, by setting 

Al = A,A e j = AJ_! + [0,e]ej := {q + te 3 : q£ A)_ u t G [0, e] }. 

Obviously, the sequence of sets Aj, j — 1,2,..., is increasing with j. The monotonicity 
property implies 

{q: $(q)<b}cA e m . 
Indeed, if < b, then for the vector q' := q — e • e we have by (2): 

$(<?') < + e • e) - e = $(<?) - e<6-e<a, 

meaning that q' G { $ < a} — A and, therefore, 

q = q + e • e G A^. 

Now, we conclude that 

{g: #(«) el} = {q: #(«) G (a, 6) } 
= {g: < 6}\{g: *(?)<a}c4\A 



Furthermore, 



For g' G R m_1 , set 



lT{q: Hl)eI}<^ m (A f m \A) 

/ m \ m 

V?'=l / j=l 



/ 1 (g') = {?i€R: (gi,g') 
By definition of set A\, this is an interval of length not bigger than e. Then we have 

fi m (At\A) = j dn m -\q') j d^ Ql ) < s(p,e). (3) 

h 

Similarly, we obtain for j = 2, . . . , m 

H m {A] \Aj_ 1 ) < s(fj,,e), 

yielding 

m 

^{q: $(g)G/}<E^ m (^\^-i)<^-^^)- - 

3 = 1 

Now, taking into account the above Remark 1, Stollmann's theorem yields immediately 
the following estimate. 
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Theorem 3 Let Ha be an LSO with random potential V(x; oS) in a finite box A C 1 d 
with Dirichlet b.c., and E(/2a) its spectrum, i.e. the collection of its eigen-values E^ A \ 
j — 1, • • • , |A|. Assume that r.v. V(x; ■) are i.i.d. with marginal distribution function Fy 
satisfying 

s(e) = sup (Fy(a + e) — Fv(a)) < oo. 

Then 

P{dist(E(ff A (w),-B) < e} < A| 2 s(e). 

3 Extension to multi-particle systems 

Results of this section have been obtained by the author and Y. Suhov [9]. 

Let N > 1 and d > 1 be two positive integers and consider a random LSO H = H(u>) 
which can be used, in the framework of tight-binding approximation, the as the Hamiltonian 
of a system of quantum particles in Z d with random external potential V and interaction 
potential U. Specifically, let xi,...,xn £ Z d be positions of quantum particles in the 
lattice Z d , and x — (xi, . . . , xn). Let {V(x; lo), x £ Z d } be a random field on 1 d describing 
the external potential acting on all particles, and U : (xi, . . . , zjv) i— > R be the interaction 
energy of the particles. In physics, U is usually to be symmetric function of its N arguments 
xi, . . . ,xn £ Z d . We will assume in this section that the system in question obeys either 
Fermi or Bose quantum statistics, so it is convenient to assume U to be symmetric. Note, 
however, that the results of this section can be extended, with natural modifications, to 
more general interactions U. Further, in [9] U is assumed to be finite-range interaction: 

supp U C {x : maxd^j — Xk \ < r)}, r < oo. 

Such an assumption is required in the proof of Anderson localization for multi-particle 
systems, however, it is irrelevant to the Wegner-Stollmann estimate we are going to discuss 
below. 

Now, let H be as follows: 

(# H/)Ge) = E ( A(J) + V{x?,uj) + U(x), 

3=1 

where is the lattice Laplacian acting on the j-th particle, i.e. 

A (i) = 1 A ® . . . ig> 1 

1 j N 

acting in Hilbert space £ 2 (Z Nd ). For any finite "box" 

A = A (1) x ... x A (Ar) c Z Nd 

one can consider the restriction, Ha{u), of H(uj) on A with Dirichlet b.c. It is easy to see 
that the potential 

JV 

w(x) = Y. v ( x r^) + u{x) 

3=1 

is no longer an i.i.d. random field on Z Nd , even if V is i.i.d. Therefore, neither Wegner's nor 
Stollmann's estimate does not apply directly. But, in fact, Stollmann's lemma does apply to 
multi-particle systems, virtually in the same way as to single-particle ones. 

Theorem 4 Assume that r.v. V(x; ■) are i.i.d. with marginal distribution function Fv 
satisfying 

s(e) — sup (Fv(a + e) — Fv(a)) < oo. 
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Then 

P{dist(E(# A (w),-E) < e} < |A| • M(A) ■ s(e), 

with 

N 

M(A) = ^ cardA (3) . 
3=1 

Proof. Fix A and consider the union of all lattice points in Z d which belong to the single- 
particle projections A"', j = 1, . . . , N: 



X(A) = \J A (J) C Z d 



3=1 

Now we can apply Stollmann's lemma to Ha by taking the index set J = X(A) and auxiliary 
probability space R J . Indeed, the random potential V(xju) := V(xi;lu) + ■ ■ ■ + V(xn',^) 
can be re-written as follows: 

V(xi;uj)-\ h V{xn;u) = ^ c(x,y)V(y;u) 

yeX(A) 

with integer coefficients c(x, y) such that 

c(x,y)>0, Yl c ^y)= N - ( 4 ) 

yeX(A) 

For example, if N = 2, one can have either V(x\,uj) + V(x2; w) with xi ^ x%, in which case 
we have 

J 1, if y = xx or y = x 2 
c{x, y) — < 

I 0, otherwise 

or V(a;i;a;) +V(a;i;aj) = 2V(a;i;w) for "diagonal" points (xi,X2), where 

/ \ f2, ify = xi 
c(x, y) = < 

I 0, otherwise 

In any case, as shows @, random potential at x S A is a linear function of one or more 
coordinates in the auxiliary space R J growing at rate > Nt > t along the principal diagonal 
{qi = 92 = ■ ■ ■ = q\j\ = t £ R}. Hence, the operators of multiplication by V(x;,uj) form a 
J-monotonic family, and, by virtue of Remark 2, the same holds for H — Ho + U + V(uj), 
just as in the single-particle case (and even "better", for N > 1 !). By Theorem 2, this 
implies immediately the estimate 

P{dist(E(tf A M,£) < e} < |A| 2 s(At,e). . 

It is not difficult to see that the same argument, with obvious notational modifications, 
applies to Fermi and Bose lattice quantum systems, i.e. to restrictions of H to the sub- 
spaces of symmetric (Bose case) or anti-symmetric (Fermi case) functions of arguments 
xi,...,x N on (Z d ) N . 

4 Extension to correlated random variables 

Now let fi m be a measure on R m with marginal distributions of order m — 1, 

m — 1 / / \ m— 1/ \ ■ -i 

AS W^i) = Mj •••,?j-l> -,9m), J = l,...,m, 
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and conditional distributions fx](qj | q'^j) on qj given all q k ,k / j. For every e > 0, define 
the following quantity: 

Ci(^i m , e) = max sup / d/i m_1 (q^) / 

J agl J J 

a 

and assume that Ci(/U, e) is finite: 

a + e 

max sup / dfj™" 1 (q'jtj) / d/i(<7i|9^j) < oo. (5) 
J a£i J J 



Remark 3 As a simple sufficient condition of finiteness o/Ci(/i,e), one can use, e.g., a 
uniform continuity (but not necessarily absolute continuity !) of the single-point condi- 
tional distributions, 



3 q'-tj «6« 



a-\-e 



max sup sup / dfi(qj\qUj) < C2(// m ,e) < oo 



or even the existence and uniform boundedness of the density p(qj\q'^j) of these conditional 
distributions: 

Bvpp(qj\q&) < C 3 (ji m ,e). 



Remark 4 In applications to localization problems, the aforementioned continuity moduli 
Ci(/i m ,e), C2(/i m ,e), C3(/i m ,e) need to decay not too slowly as e — > 0. A power decay of 
order 0(e 13 ) with (5 > is certainly sufficient, but it can be essentially relaxed. For example, 
it suffices to have an upper bound of the form 

Ci (/x m ,e~ L ' 3 ) < Const -IT 3 , 

uniformly for all sufficiently large L > with some (arbitrarily small) (3 > and with B > 
which should sufficiently big, depending on the specific spectral problem. 

Using notations of the previous section, one can formulate the following generalization of 
Stollmann's lemma. 

Lemma 1 Let <£> : R J — > R, E J = E m , be a J-monotonic function and fi m a probability 
measure o»l m = R J with Ci(/i m , e) < oo. Then for any interval I C R of length \I\ = e > 0, 
we have 

lT{q: <!>(<?) €/}<m-Ci(M,e). 

Proof. We proceed as in the proof of Stollmann's lemma and introduce in R m the sets 
A — { q : $(q) < a } and Aj, j = 0, . . . , m. Here, again, we have 

{q: $(q) G / } =C A c m \A 

and 



/**"{«: $( ? )G/}<^M m {A)\A)_ X ) 

3 = 1 
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For q'-t! G R m ~\ we set 

h(q'^ 1 ) = {qi GR: (qi, q'^) G A\ \ A }. 
Furthermore, we come to the following upper bound which generalizes (3): 

M m (^ \ A) = | d M m ~V) | rfM(aik') < CrOu.e). (6) 
Similarly, we obtain for j = 2, . . . ,m 

M m (A5\^_o <Ci (/.,£), 

yielding 

m 

5 Application to Gaussian random fields 

Let V(£, w),a; G Z d , d > 1, be a regular stationary Gaussian field of zero mean on the lattice 
Z d . The regularity implies that the field V(-,ui) is non-deterministic, i.e. the conditional 
probability distribution of 1^(0,-) given {V(y),y 7^ 0} is Gaussian with strictly positive 
variance. In other terms, the r.v. V(0, •), considered as a vector in the Hilbert space H. VZ d 
generated by linear combinations of all V(x, •), x G Z d , with the scalar product 

K,J/)=E[£!,], 

does not belong to the subspace T~Lv i. d \{o}'- 

\\v(o,-) -n HvzdKW v(o,-)f = a% > 0, 

where 

Un v ^ = ^\V(x,-),xeI. d \{Q}]. 
Furthermore, for any subset A C Z d \ {0}, 

\\v(o,-)~n nv . A v(o,-)\\ 2 >tl 

since 

Tiv.A C Ti.v,z d \{0}- 

Therefore, the conditional variance of V(0, •) given any non-zero number of values of V 
outside x — is bounded from below by a 2 . Respectively, the conditional probability density 
of V(0, ■), for any such nontrivial condition is uniformly bounded by (27ro"o) -1 ^ 2 < 00. Now 
a direct application of Lemma [1] leads to the following statement. 

Theorem 5 Let A C 1 d be a finite subset of the lattice, and A' C 1 d \ A any subset disjoint 
with A (A' may be empty). Consider a family of LSO Ha(u) with Gaussian random potential 
V{ijj) in A, with Dirichlet b.c. on dA. Then for any interval I C R of length e > 0, we have 

P { E(# A ) n I # I V(y, •), y G A' } < G(V) |A| 2 e, 
where the constant C(V) < 00 whenever the Gaussian field V is non-deterministic. 
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6 Application to Gibbs fields with continuous spin 



Apart from Gaussian fields, there exist several classes of random lattice fields for which the 
hypothesis of Lemma[T]can be easily verified. For example, conditional distributions of Gibbs 
fields are given explicitly in terms of their respective interaction potentials. Specifically, 
consider a lattice Gibbs field s(x,u) with bounded continuous spin, 

s : Q x Z d -> S = [a, b] C R 

generated by a short-range, bounded, two-body interaction potential u(-, •). The spin space 
is assumed to be equipped with the Lebesgue measure ds. In other words, consider the 
formal Hamiltonian 

H(s) = h ( x ) + "l*-y|0O),w(i/)), 

i6Z d xgZ d \v-x\<R 

where h : S — > R is the self-energy of a given spin. The interaction potentials u\ x _ y \ (s(x),s(y)) 
vanish for \x — y\ > R and are uniformly bounded: 

max sup |u;(s,i)| < oo. 

1 < R a,tes 

Then for any lattice point x and any configuration s' = s'u. x of spins outside {x}, the single- 
site conditional distribution of s(x) given the external configuration s' admits a bounded 
density 

e -0U(s*\s'\ e -PU(.s*\8'\ 

P( s * I = H(As /) = j s e-0u { t\s' )dt 

with 

U{s x \s) := ^2 u\ x _y\(s x ,s' y ) 

y: \y — x\ <R 

satisfying the upper bound 

\U(s x \s')\ < (2R + l) d sup \ui(s,t)\ < oo. 
s,tes 

A similar property is valid for sufficiently rapidly decaying long-range interaction potentials, 
for example, under the condition 

sup \u M (s,t)\< °™l s , 5>0. (7) 

as well as for more general, but still uniformly summable many-body interactions. Here is 
one possible Wegner-Stollmann-type result concerning such random potentials. 

Theorem 6 Let A C Z d be a finite subset of the lattice, A'c2 i \A any subset disjoint with 
A (A may be empty), and let s(x, u>) be a Gibbs field in A with continuous spins s £ S — [a, b] 
generated by a two-body interaction potential Ui(s,t) satisfying condition Q), with any b.c. 
on Z d \ A. Consider a LSO Ha with random potential V(x,u>) = s(x,lo). Then for any 
interval I C R of length e > 0, we have 

P { E(H A ) I # | V{y, •), V 6 A' } < C(V) \A\ 2 e, C(V) < oo. 

In the case of unbounded spins and/or interaction potentials, the uniform boundedness of 
conditional single-spin distributions does not necessarily hold, since the energy of interaction 
of a given spin s(0) with the external configuration s' may be arbitrarily large (depending on 
a particular form of interaction) and even infinite, if s'(y) — > oo too fast. In such situations, 
our general condition (J5J may still apply, provided that rapidly growing configurations s' 
have sufficiently small probability, so that the outer integral in the r.h.s. of © converges. 
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7 Conclusion 



Wegner-Stollmann-type estimate of the density of states in finite volumes is a key ingredient 
of the MSA of spectra of random Schrodinger (and some other) operators. The proposed 
simple extension of Stollmann's lemma shows that a very general assumption on correlated 
random fields generating potential rules out an abnormal accumulation of eigen-values in 
finite volumes. This extension applies also to multi-particle systems [10] . 
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